Spectral signatures of Holstein polarons 

HOLGER Fehske and Andreas Alvermann 

Institut fiir Physik, Ernst-Moritz-Arndt-Universitdt Gretfswald, 17487 Greifswald, Germany 
Martin Hohenadler 

Institut fiir Theoretische Physik ~ Computational Physics, TU Graz, 8010 Graz, Austria 
Gerhard Wellein 

Regionales Rechenzentrum Erlangen, Universitat Erlangen, 91058 Eriangen, Germany 



1. Fundamentals 

I'l. Self-trapping phenomenon. - Electrons or holes delocalised in a perfect rigid 
lattice can be "trapped" in a potential well produced by displacements of atoms if 
the particle-lattice interaction is sufficiently strong [1]. Such trapping is energetically 
favoured over wide-band Bloch states if the carrier's binding energy exceeds the strain 
energy required to produce the trap. Since the potential itself depends on the carrier 
state, this highly non-linear process is called "self-trapping" (Fig. 1). A self-trapped 
state is referred to as "large" if it extends over multiple lattice sites. Alternatively, for 
a quasi-particle (QP) with extremely large effective mass to* - practically confined to a 
single site - the state is designated as "small" . Self-trapping does not imply a breaking 
of translational invariance, i.e., in a crystal these eigenstates are still itinerant allowing, 
in principle, for coherent transport with an extremely small bandwidth. 

Introducing the concept of polarons into physics, the possibility of electron self- 
trapping was pointed out by Landau as early as 1933 [2]. Self-trapped polarons consisting 
of electrons accompanied by phonon clouds can be found, e.g., in alkali(ne) metal (earth) 
halides, II-IV- and group- IV semiconductors, and organic molecular crystals [3]. With the 
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Fig. 1. - Polaron formation: The 
standard phase transition concept 

fails to describe this highly non- 
linear feedback phenomenon. 



observation of polaronic effects in high- Tc cuprates [4] and colossal magneto-resistance 
manganites [5], research on polarons has attracted renewed attention (see also [6]). 

1'2. Holstein model. - Depending on the relative importance of long- and short-range 
electron-lattice interaction, simplified models of the Frohlich [7] or Holstein [8] type have 
been widely used to analyse polaronic effects in solids with displaceable atoms. The 
spinless Holstein Hamiltonian reads 

(1) ^ = ~* S '''''^ ~ V^Yl^bl + b,)n, +u;oYl ^Ih • 

Here c'^' and b^^^ are the annihilation [creation] operators of a (spinless) fermion and a 
phonon at Wannicr site i, respectively, and = c[c^. In (1), the following idealisations 
of real electron-phonon (EP) systems are made: (i) the electron transfer t is restricted to 
nearest-neighbour pa,irs {ij); (ii) the charge carriers are locally coupled to a dispersionless 
optical phonon mode [e^ measures the polaron binding energy and uq is the bare phonon 
frequency (h = 1)]; (iii) the phonons are treated within the harmonic approximation. 

Nevertheless, as yet, none of the various analytical treatments, based on weak- and 
strong-coupling adiabatic and anti-adiabatic perturbation expansions [9] , are suited to in- 
vestigate the physically most interesting polaron transition region. In the latter, the char- 
acteristic electronic and phononic energy scales are not well separated and non-adiabatic 
effects become increasingly important, implying a breakdown of the standard Migdal 
approximation. Quasi-approximation-free numerical methods like quantum Monte Carlo 
(QMC) [10-12], exact diagonalisation (ED) [13], or the density matrix renormalisation 
group (DMRG) [14] can, in principle, close the gap between the weak- and strong-EP- 
coupling limits, therefore representing the currently most reliable tools to study polarons 
close to the cross-over regime (see also [15]). 

1"3. Ground-state properties. - Previous numerical work mainly focused on the single- 
polaron problem. There are two control parameters common in use, the first being the 
adiabaticity ratio u)o/t. In the adiabatic limit u>o/t 0, the motion of the particle 
is affected by quasi-static lattice deformations, whereas in the opposite anti-adiabatic 
limit LOo/t — » 00, the lattice deformation adjusts instantaneously to the carrier position. 
The second parameter is the EP coupling strength A = ep/2Dt, the ratio of the polaron 
binding energy of an electron confined to a single site and the free electron half-bandwidth 
in D dimensions. 

Figure 2 summarises selected ground-state characteristics of small-polaron formation, 
showing results obtained by (variational) ED [16,17] and DMRG [14]. The transition 
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Fig. 2. - Effective mass m* /mo (data taken from [16]), (normalised) electron-phonon correlation 
function xo,j = ('T'o(&o+j +&o+j))/25('T'o)5 and schematic "phase diagram" of the Holstein model. 
Depending on the adiabaticity of the system, the cross-over regime is determined by the more 
stringent of the two conditions A ~ 1 and gr ~ 1. Thus starting from "light" {uJo/t < 1) 
or "heavy" {uio/t > 1) electrons it is possible to understand the formation of small adiabatic 
"Holstein" or anti-adiabatic "Lang-Firsov" polarons as two limiting cases of a general picture. 



to a polaron with large effective mass takes place at A ~ 1, and is much sharper in 
higher dimensions [16]. Adiabatic theory predicts an energy barrier separating quasi-free 
(infinite radius) polarons from small-sized lattice polarons in D > 1, but not in D = 1. 
However, the (variational) ED data show that for ojq > 0, the cross-over is continuous in 
any dimension, i.e., the Holstein model does not exhibit a true phase transition [16, 17], 
in accordance with the theorem of Gerlach and Lowen [18]. In the anti-adiabatic regime, 
a further parameter ratio, — Ep/ujQ, turns out to be crucial. Small-polaron formation 
takes place if both A > 1 and > 1. Note that the mass enhancement notably deviates 
from the strong-coupling (perturbation theory) result exp((7^) at intermediate phonon 
frequencies [16]. The cross-over from an only weakly phonon-dressed electron to a small 
polaron also shows up very clearly in the correlation between electron- and phonon- 
displacements (not be equated with the polaron radius): At the critical coupling, there is 
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a strong enhancement of the on-site correlations [14, 17]. Furthermore, Xo,j/Xo,o decays 
more rapidly as D increases, i.e., the surrounding phonons are localised closer to the 
electron in higher dimensions [16]. 

The ground-state properties discussed so far constitute only one aspect of polaron 
physics. Of equal importance is the dynamical response of a polaronic system to external 
perturbations. For a moving (small) polaron there is a perpetual exchange of momentum 
between the electron and the deformation field [19]. Photoemission spectroscopy and 
inelastic neutron scattering should therefore be able to detect the strong interrelation of 
electron and phonon degrees of freedom. The strong mixing of purely electronic states 
with lattice vibrational excitations should be visible in the optical response as well. In 
the following sections, we therefore present exact numerical results for clcctron/phonon 
spectral functions and the optical conductivity in the framework of the Holstein model. A 
very efficient Chebyshev-expansion-based algorithm for the calculation of such dynamical 
correlation functions is outlined in [15] (for more details see [20]). 

2. — Photoemission spectra 

Examining the dynamical properties of polarons, it is of particular interest whether a 
QP-like excitation exists in the spectrum. This is probed by direct (inverse) photoemis- 
sion, where a bare electron is removed (added) from (to) the many-particle system. The 
intensities (transition amplitudes) of these processes are determined by the imaginary 
part of the retarded one-particle Green's functions, 

(2) G±(fc,c^) = ((cJ;c±))^= hm (Vo|cn'^ + i^Tif]-^c±l^o), 

r;— »0+ 

i.e., by the wave- vector resolved spectral functions 

(3) A^{k,u,) = -ilmG±(fc,a;) = ^ l(V^lc±l^o)l' ^[^ T {E^ - E,)] 

m 

with = cl and = Cfe- These functions test both the excitation energies E^ — Eq and 
the overlap of the ground state \ipo) with the exact eigenstates \tp^) of a {Ng ± l)-particle 
system. Hence, G+ {k, to) [G~ {k,Lo)] describes the propagation of an additional electron [a 
hole] with momentum k [—k] and energy w. The electron spectral function of the single- 
particle Holstein model corresponds to = 0, i.e., A~{k,uj) = 0. A{k,uj) = A~^{k,uj) + 
A~{k,uj) can be determined, e.g., by cluster perturbation theory (CPT) [15,21]: We first 
calculate the Green's function G^j (to) of a A^c-site cluster with open boundary conditions 
for i,j = 1, . . . -Nc, and then recover the infinite lattice by pasting identical copies of 
this cluster along the edges, treating the inter-cluster hopping in first-order perturbation 
theory. 

Figure 3 shows that at weak coupling (left panel), the electronic spectrum is nearly un- 
affected for energies below the phonon emission threshold. Hence, for the case considered 
here with ojq lying inside the bare electron bandwidth ADt, the renormalised dispersion 
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Fig. 3. - Spectral function of the ID liolstcin polaron calculated within CPT in the weak (a) 
and strong (b) non-adiabatic EP coupling regime. CPT is based on ED of a finite cluster with 
Nc sites, M = 7 (A = 0.25) and M = 25 (A = 2) phonon quanta. 



E{k) follows the tight-binding cosine dispersion (lowered 3c e,,) up to some kx, where the 
dispersionless phonon intersects the bare electron band. For k > kx, electron and phonon 
states "hybridise" , and repel each other, leading to the well-known band-flattening phe- 
nomenon [22]. The high-energy incoherent part of the spectrum is broadened oc Sp, with 
the fc-dependent maximum again following the bare cosine dispersion. 

To substantiate this interpretation, we compute the QP weight and mean phonon 
number, 

(4) Zik) = \{^+\cl\%)\' and N^ik) = \{bl + b,)\%) , 

for each /c-sector (see Fig. 4). The fc-dependent Z-factor can be taken as a measure of 
the "electronic contribution" to the QP. For weak EP coupling, we have Z{k <C fc^) ^ 1 
[Z{k > fcx) <C 1], reflecting the electronic [phononic] character of the states at the band 




Fig. 4. - QP residue Z(k) and mean phonon 
number Np{k) (inset) obtained by ED at the 
allowed k vectors of a ID 14-site lattice with 
periodic boundary conditions. Note that 
roughly gives the mean phonon number in the 
strong-coupling case. 
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centre [band edge], which are basically zero-phonon [one-phonon] states (see Fig. 4, inset). 
Increasing the EP coupling, a strong mixing of electron and phonon degrees of freedom 
takes place, whereby - forming a small polaron - both quantum objects cease to exist 
independently. As expected, this leads to a significant suppression of the (electronic) 
QP residue Z{k) for all k. Now the states \tpo)j li'k) multi-phonon states, and the 
polaronic QP is heavy because it has to drag with it a large number of phonons in its 
phonon cloud (cf. Fig. 4, inset). 

The inverse photoemission spectrum in the strong- coupling case is shown in the right 
panel of Fig. 3. First, we observe all signatures of the famous polaronic band-collapse, 
where a well-scparatcd, narrow (i.e., strongly renormalised) , coherent QP band is formed 
at a; ~ —Bp. If wc had calculated the polaronic instead of the electronic spectral func- 
tion (3), nearly all spectral weight would reside in the coherent part, i.e., in the small- 
polaron band [23]. In contrast, Z{k), defined by (4), is extremely small and approaches 
the strong-coupling result Z =- cxp(— g^) for A , 1. Note that the inverse effective 

mass m* /rriQ and Z{k) differ if the self-energy is strongly A;-dependent. This discrepancy 
has its maximum in the intermediate-coupling regime for ID systems, but vanishes in 
the limit A ^ oo and, in any case, for D = oo [16]. The incoherent part of the spectnun 
is split into several sub-bands separated in energy by ujq, corresponding to excitations of 
an electron and one or more phonons (Fig. 3). 

Finally, let us emphasise that for all couplings, the lowest-lying band in A{k, w) almost 
perfectly coincides with the coherent polaron band-structure (solid line in Fig. 3) obtained 
by variational ED [16]. This method has been shown to give very accurate results for 
the infinite system, also underlining the high precision of our CPT approach. 

3. — Phonon spectra 

Next we show that the phonon spectra provide additional useful information concern- 
ing the polaron dynamics. For this purpose, we calculate the T = spectral function 
B{q,uj), which is related to the phonon Green's function D{q,uj) for w > by 

(5) B{q,w) = --lm.D{q,w) , D{q,u)) = 2u!o{{xq; x_q))^ , 

TT 

where Xq = N~2 ^jXje~^^'' and Xj = (6j -|- bj)/\/2ujQ. Again we employ a cluster 
approach, determining first the cluster phonon Green's fimction D|^ (tj), and afterwards, 
as in CPT, constructing the energy and momentum dependent Green's function of the 
infinite lattice from D{q,u)) = X]fj=i ^iji'^) e~'*^'~-'^. Since for the Holstein model 
the bare phonon Green's function £'°(w) = 2ujo/{uj^ — u)q) is fc-independent, this cluster 
expansion is identical to replacing the full real-space Green's function D^- by D?-. 

Figures 5(a) and (b) show the evolution of the phonon spectrum with increasing 
EP interaction in the adiabatic case (wo/f = 0.4). For A = 0.5, there is a dispersive 
low-energy absorption reflecting the polaron band dispersion, as demonstrated by the 
comparison with variational-ED data for E{k). If the weakly renormalised electron band 
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Fig. 5. - Phonon spectral function for the 
ID single-electron Holstein model. The bare 
phonon frequency is indicated by a dashed 
vertical line. The solid line corresponds to 
the polaron band dispersion E{k) — E{0), cal- 
culated using the method presented in [16]. 
Shown are exemplary results for the adiabatic 
weak and strong-EP-coupling [(a) and (b)] and 
anti-adiabatic intermediate-coupling [panel (c)] 
regimes. 



intersects the (dispersionless) bare phonon excitation at some qy, level repulsion occurs 
and we observe two adjacent absorption features. At larger EP couplings, the polaron 
band separates completely from the bare phonon signature until, in the extreme strong- 
coupling limit, two almost flat absorption bands emerge, corresponding to the lowest 
(small) polaron band and the first excited band, separated by a one-phonon excitation. 
Both signals are strong because of the large "phonon content" in these states. 

As stated above, in the anti-adiabatic regime, the small-polaron cross-over is deter- 
mined by the ratio g^, and occurs at about — 1. For 5^ <C 1, nearly the whole spectral 
weight resides in the bare phonon peak (if ujq > 4Z3t). The phonon spectrum near the 
transition point is shown in Fig. 5(c) for ujo/t — 4. We detect a clear signature of the 
small-polaron band with renormalised width W 1.5i - a factor of about ten larger 
than in the adiabatic case {X = 1, u!o/t — 0.4; note that the excitation energy in Fig. 5 
is plotted in units of luq). The dispersionless excitation at ujq is obtained by adding one 
phonon with momentum q to the fc = ground state. Above this pronounced peak, 
we find an "image" of the lowest polaron band - shifted by ujq - with extremely small 
spectral weight, hardly resolved in Fig. 5. 
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4. Optical response 

We apply the ED-KPM scheme outhned in [15, 20, 24] to calculate the optical absorp- 
tion of the ID Holstein polaron. The results for the (regular) real part of the conductivity. 



(6) 



E,„>Eo 



(here j — ~^e't^i{clc.^_^_i — c|^_,_]^Cj) is the current operator), and possible deviations from 
established polaron theory are important for relating theory to experiment. The standard 
description of small polaron transport [25] yields for the T = ac conductivity 



(7) 



RetT(w) 



(To 1 



■ exp 



4£pWo 



For sufficiently strong coupling this formula predicts a weakly asymmetric Gaussian 
absorption peak centred at twice the polaron binding energy. 

Figure 6 shows Rei7(w) when polaron formation sets in (a), and above the transition 
point (b). For A = 2 and coo/t = 0.4, i.e., at rather large EP coupling but not in the 
extreme small-polaron limit, we find a pronounced maximum in the low-temperature 
optical response, which, however, is located somewhat below 2ep = 2g^ijjQ, the value for 
small polarons at T = 0. At the same time, the line-shape is more asymmetric than in 
small-polaron theory, with a weaker decay at the high-energy side, fitting even better to 
experiments on standard polaronic materials such as Ti02 [26]. At smaller couplings, 
significant deviations from a Gaussian- like absorption are found, i.e., polaron motion is 
not adequately described as hopping of a self-trapped carrier almost localised on a single 
lattice site. 




0.5 1 
CO/ 28 



0.5 1 
CO/ 28 



Fig. 6. - Optical conductivity (in units of Tre^t^) of the ID Holstein polaron at T = compared 
to the analytical small-polaron result (dashed line). ED data are for a system with N = 6 sites 
and 45 phonons; tjo is determined to give the same integrated spectral weight as Kea(ijj > 0). 
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Fig. 7. - One-electron spectral function A{k,Lj — fi) of the ID spinless Holstein model from 
QMC [11] for low (n = 0.1) and high (n = 0.4) carrier densities at weak [(a),(b); A = 0.1] and 
strong [(c), (d); A = 2] EP coupling, uio/t = 0.4, A'^ = 32, and t/ksT = 8 (inverse temperature). 
Note that QMC cannot resolve features with very small spectral weight. 



5. Many-polaron problem 

Finally, we address the important issue how the character of the QPs of the system 
changes with carrier density. While for very strong EP coupling no significant changes 
are expected due to the existence of rather independent small (self-trapped) polarons 
with negligible residual interaction, a density-driven cross-over from a state with large 
polarons to a metal with weakly dressed electrons may occur in the intermediate-coupling 
regime [11]. This problem has recently been investigated experimentally by optical mea- 
surements on La2/3(Sr/Ca)i/3Mn03 films [27]. 

In the Holstein model, the above-mentioned density-driven transition from large po- 
larons to weakly EP-dressed electrons is expected to be possible only in ID, where large 
polarons exist at weak and intermediate coupling. The situation is different for Frohlich- 
type models [7, 28, 29] with long-range EP interaction, in which large-polaron states exist 
even for strong coupling and in D > 1. 

To set the stage, Fig. 7 shows the evolution of the one-electron spectral function 
A{k,uj — /i) with increasing electron density n in the weak- and strong- coupling limiting 
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Fig. 8. - Single particle spectral functions A~{k,u!) (red dashed lines) and A'^{k,u}) (black solid 
lines) for different band fillings n at uio/t = 0.4 and A = 1 (CPT results with Nc — 10). Blue 
crosses track the small-polaron band obtained by ED. 



cases. In the former [Fig. 7(a), (b)], the spectra bear a close resemblance to the free- 
electron case for all n, i.e., there is a strongly dispersive band running from —2t to 
2t, which can be attributed to weakly dressed electrons with an effective mass close 
to the non-interacting value. As expected, the height (width) of the peaks increases 
(decreases) significantly in the vicinity of the Fermi momentum, which is determined by 
the crossing of the band with the chemical potential ^. In the opposite strong-coupling 
limit [Fig. 7(c),(d)], the spectrum exhibits an almost dispersionless coherent polaron band 
at fi (QMC, due to the use of maximum entropy, has problems resolving such extremely 
weak signatures). Besides, there are two incoherent features located above and below 
the chemical potential, broadened oc e^, reflecting phonon-mediated transitions to high- 
energy electron states. At n = 0.4, the photoemission spectrum for k < tt/2 becomes 
almost symmetric to the inverse photoemission spectrum for k > tt/2 and already reveals 
the gapped structure expected at n — 0.5 due to the Peierls transition (charge-density- 
wave formation). The most important point, however, is the clear separation of the 
coherent band from the incoherent parts even at large n. This indicates that small 
polarons are well-defined QPs in the strong-coupling regime, even at high carrier density. 

Figure 8 displays the inverse photoemission [A+(fc,a;)] and photoemission spectra 
[A~{k,u!)] at intermediate EP coupling strength, determined by CPT at T = 0. 

At n = 0.1 we can also identify a (coherent) polaron band (blue crosses; cf. also the 
ED data in [11]). The Fermi energy E-p = /i(T 0) is located within this signature. 
The large-polaron band has small electronic spectral weight especially away from E-p and 
fiattens as an effect of the EP coupling at large k (see Sec. 2). Below this band, there 
exist equally spaced phonon satellites, reflecting the Poisson distribution of phonons 
in the ground state. Above E-p there is a broad dispersive incoherent feature whose 
maximum follows closely the dispersion relation of free particles. 
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As n increases, a well-separated coherent polaron band can no longer be identified. 
At about n ~ 0.3 the deformation clouds of the (large) polarons start to overlap leading 
to a mutual (dynamical) interaction between the particles. Increasing the carrier den- 
sity further, the polaronic QPs dissociate, stripping their phonon cloud. Now diffusive 
scattering of electrons and phonons seems to be the dominant interaction mechanism. 
As a result both, the phonon peaks in A~{k,u}) and the incoherent part of A'^{k,u}) are 
washed out, the spectra broaden and ultimately merge into a single wide band [11,30]. 
As can be seen from Fig. 8 (b), the incoherent excitations lie now arbitrarily close to the 
Fermi level. These are significant differences to the nearly free electron and small polaron 
spectra shown for the same carrier density in Figs. 7 (b) and (d), respectively. Obviously 
the low-energy physics of the system can no longer be described by small-polaron theory. 

6. — Summary and open problems 

In this contribution, we have reviewed ground-state- and, most notably, spectral prop- 
erties of Holstein polarons by means of quasi-exact numerical methods such as (varia- 
tional) Lanczos diagonalisation, a kernel polynomial expansion technique, cluster per- 
turbation theory, and quantum Monte Carlo. Our numerical approaches yield unbi- 
ased results in all parameter regimes, but are of particular value in the non-adiabatic 
intermediate-coupling regime, where perturbation theories and other analytical tech- 
niques fail. The data presented for the (inverse) photoemission, phonon- and optical 
spectra show that electron and phonon excitations become intimately, dynamically con- 
nected in the process of polaron formation. 

Although wc have now achieved a rather complete picture of the single (Holstein) 
polaron problem (perhaps dispersive phonons, longer-ranged EP interaction, finite tem- 
perature and disorder effects deserve closer attention), the situation is discontenting in 
the case of a finite carrier density. Here electron-phonon coupling competes with some- 
times strong electronic correlations as in, e.g., unconventional ID metals, quasi-lD MX 
chains, quasi-2D high- Tc superconductors, 3D charge-ordered nickelates, or bulk colos- 
sal magneto-resistance manganites. The corresponding microscopic models contain (ex- 
tended) Hubbard, Heisenberg or double-exchange terms, and maybe also a coupling to 
orbital degrees of freedom, so that they can hardly be solved even numerically with the 
same precision as the Holstein model. Consequently, the investigation of these materials 
and models will definitely be a great challenge for solid-state theory in the near future. 

* * * 

We would like to thank J. Bonca, D. Ihle, J. Loos, S. A. Trugman, G. Schubert and 
A. Weifie for valuable discussions. 
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